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We emphasize the importance of the position dependence of the diffusion coefficient D{r) in 
the self-consistent theory of localization and argue that the scaling law T oc Ini/L^ obtained by 
Cheung and Zhang [Phys. Rev. B 72, 235102 (2005)] for the average transmission coefficient T 
of a disordered slab of thickness L at the localization transition is an artifact of replacing D{r) by 
its harmonic mean. The correct scaling T oc l/L^ is obtained by properly treating the position 
dependence of D{y). 

PACS numbers: 42.25.Dd 



In a recent paper Cheung and Zhang (CZ) apply 
the self-consistent (SC) theory of localization to study 
the transmission of waves through a slab of disordered 
medium at the Anderson localization transition. The SC 
theory is a powerful tool to deal with the phenomenon of 
Anderson localization, but its application to disordered 
media of finite size requires some care. In the original 
papers by VoUhardt and Wolfle the size L of dis- 
ordered sample was acknowledged using a lower cut-off 
in the integration over momentum. Despite the obvious 
crudeness of this approach, it was sufficient to recover 
the main results of the scaling theory of localization [3] 
and added a great physical insight into the phenomenon 
of disorder-induced localization. Later on. Van Tigge- 
len et al. U argued that in a medium of finite size the 
SC theory naturally leads to a position dependence of 
the diffusion coefficient D{y). This adapted SC theory 
was successfully applied to study coherent backscatter- 
ing [4, and dynamics jSJ [5] of localized waves. Micro- 
scopic justifications for position dependence of D have 
been recently presented based on the diagrammatic [7 
and field-theoretic [8j calculations. 

CZ propose a way of overcoming technical difficulties 
caused by the position dependence of D{t) (see also 
[H]). They average the equation for 1/Z)(r), their Eq. (1), 
over the sample volume, thus replacing D{y) by its har- 
monic mean D. In this Comment we argue that although 
such an approach can be justified in the weak localization 
regime [S], it is not adequate at the mobility edge and in 
the Anderson localization regime. In particular, our cal- 
culations that properly treat the position dependence of 
D(r), do not confirm the scaling law T cx In L/L^ found 
by CZ for the transmission coefficient T of a disordered 
slab of thickness L at the mobility edge. Instead, we 
find T cx l/L^ in agreement with the scaling theory of 
localization 

To study the scaling of the average transmission coeffi- 
cient T with the thickness L of disordered slab, we solve 
the two equations of SC theory — Eqs. ( 1) and (2) of Ref. 
[6] with n = (stationary regime) and M — \ (mobility 
edge) |10j — numerically. We use the same boundary 
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FIG. 1: Average transmission coefficient T of a disordered 
slab of thickness L at the Anderson localization transition. 
Circles were obtained from the self-consistent theory of local- 
ization with a position-dependent diffusion coefficient D[z) by 
numerical solution . The solid red line is a fit to the numeri- 
cal results using Eq. ^ with D{0)/Db = 0.82 and Zc = 4.2^. 
The dotted blue and dashed green lines are fits to numeri- 
cal data for L/£ > 10^ using T oc {£/Lfln{L/£) and T oc 
{e/Lf ln(L/a£), respectively. We obtain a ~ 7.37 x 10"^^ in 
the latter case. 



conditions and the same method of numerical solution as 
in Ref. [6 and vary the thickness of the slab L from 10^^ 
to 8 X lO^e. Here k is the wave number of the wave and 
£ is the mean free path due to disorder. Our results are 
presented in Fig. [l] by circles. The red solid line in Fig. 
[1] shows 
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that we obtained by assuming D{z) — D{0)/{1 + z/zc 
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with z = min(z, L — z) as suggested by Van Tiggelen et 
al. (4j . Here I?b is the diffusion coefficient in the absence 
of macroscopic interferences (i.e. in the limit oi k£ ^ 1). 
D(0) / Db was determined directly from the numerical re- 
sults at a sufficiently large L = fO^^, whereas Zc was a 
free fit parameter. We used zq = f ^, corresponding to 
no internal reflections at the sample boundaries. Devia- 
tions of the fit from the numerical results do not exceed 
3% in the whole range of considered L's, which supports 
the validity of Eq. ([T|) and its underlying model for D{z). 
The inaccuracy of the latter model in the middle of the 
slab cause deviations at small L < 10^^, whereas devi- 
ations at large L > 4 x 10'^^ are mostly due to the ex- 
tremely slow convergence of our computational algorithm 
for thick slabs and would, most likely, disappear if more 
computer time were available. We note that T x (L/i)"^ 
grows with ln{L/£) for L < 10'^ £, but then saturates at a 
constant level for larger L, suggesting T cx (i/L)'^ in the 
limit of large L. 

Neither the ensemble of numerical results of Fig. [T] 
nor its small- or large-L parts can be fit by T = const x 
{£ / L)"^ ln{L / £) proposed by CZ. This is easy to see from 
Fig. [T] where we show a fit of the above equation to our 
numerical data for L/£ > 10^ (dotted blue straight line). 
It is clear that the fast growth of T x (L/f)^ with ln(L/^) 
predicted by CZ is not supported by our numerical cal- 
culations: the numerical results only show an increase of 



20% in the range of L/£ = 100-8000 and 4% in the range 
L/£ = 1000-8000, whereas the result of CZ increases by 
100% and 30%, respectively. For large L > 10^^, a rea- 
sonable fit can be achieved by T c>c {£/L)'^ ln{L/a£). The 
result of CZ would correspond to a ^ 1 , whereas a fit to 
the numerical data yields a ^ 10^^"* <C 1. This value is 
unphysically small and implies existence of length scales 
that are 24 orders of magnitude shorter than the mean 
free path £. We therefore conclude that our numerical 
results exclude the possibility of logarithmic scaling of 
T X {L/£)^ with L/£. Appearance of this scaling in Ref. 
[T] should then be an artifact of replacing D{r) by its 
harmonic mean. 

In conclusion, we have shown the importance of prop- 
erly treating the position dependence of the diffusion co- 
efficient D{r) in the SC theory of localization. In par- 
ticular, replacing D(r) by its harmonic mean leads to an 
incorrect scaling law for the transmission coefficient T 
with the thickness L of disordered slab at the mobility 
edge. The correct scaling law T cx 1/L^ is obtained by 
solving SC equations with a position dependent D{r). 

The computations presented in this paper were per- 
formed on the cluster HealthPhy (CIMENT, Grenoble). 
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ANR (project No. 06-BLAN-0096 CAROL) and the 
French Ministry of Education and Research. 
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